Anomalous Josephson Current in Junctions with Spin-Polarizing Quantum Point Contacts 
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We consider a ballistic Josephson junction with a quantum point contact in a two-dimensional electron gas 
with Rashba spin-orbit coupling. The point contact acts as a spin filter when embedded in a circuit with normal 
electrodes. We show that with an in-plane external magnetic field an anomalous supercurrent appears even 
for zero phase difference between the superconducting electrodes. In addition, the external field induces large 
critical current asymmetries between the two flow directions, leading to supercurrent rectifying effects. 

PACS numbers: 74.45.+C, 71.70.Ej, 72.25.Dc, 74.50.+r 



Josephson junctions (JJ) are the basic building blocks for 
superconducting electronics with applications that range from 
SQUID magnetometers to possible quantum computing de- 
vices. In superconductor-normal metal-superconductor (S- 
N-S) junctions the supercurrent flow is due to the Andreev 
states — a coherent superposition of electron and holes states. 
These states depend on the electronic structure of the normal 
material and on the properties of the S-N interface HI 0,(1,13]. 
Modern technologies based on two dimensional electron gases 
(2DEGs) 0] or nanowires J3] allow for a precise control 
of such electronic properties, and thus of the JJ characteris- 
tics. Moreover, spin-orbit (SO) effects offer new alternatives 
to control the spin and charge transport Jf! |9[] • 

Superconducting rectifiers are among the new devices pro- 
posed and studied during the last few years. Most of these pro- 
posals are based on the dynamics of vortices EE!- Here we 
show that in systems with SO-coupling rectifying properties 
can be obtained by controlling the spin of the Andreev states. 
To this end we consider a ballistic JJ with a quantum point 
contact (QPC) in a 2DEG with SO interaction. The QPC can 
be tuned to control the number of transmitting channels and 
thus the critical current of the junction H, 0, fi E3- B • On 
the other hand, the QPC with SO coupling may act as a spin 
filter producing spin-polarized currents when embedded in a 
circuit with normal leads Iffl U3. flol . The normal current 
also generates an in-plane magnetization — perpendicular to 
the current — as well as out-of -plane spin-Hall textures lfl7ll . 
Both effects are maximized at the core of the QPC 0]. As 
the SO-coupling preserves time-reversal symmetry (TRS), we 
expect that these peculiarities of the transmitting channels do 
not harm the Josephson effect when the leads become super- 
conducting. However, the Josephson current itself breaks the 
TRS and, as we show below, it reveals striking effects of the 
SO-coupling. For example, the supercurrent generates spin 
polarization in the 2DEG 1 19] and the QPC in a similar way 
normal current does 117, I18I1 . This is due to the distinctive 
spin texture of each Andreev state that contributes to the local 
magnetization in a supercurrent-carrying state. 

More striking effects take place if an external in-plane mag- 
netic field is applied. Its effect on the supercurrent character- 
istics depends on the nature of the junction. In the absence 



of SO-coupling, the Zeeman field may generate 7r-junctions 
resembling the case of S-ferromagnet-S junctions 111 12011 . For 
systems with SO-coupling, the existing theories include the 
description of perfectly contacted 2DEG junctions [21], wide 
junctions |22l|. I D-conductors [23] and junctions with quan- 
tum dots B24 12511 . In our case, the QPC, the internal SO field 
and the external Zeeman field conspire to reveal novel effects. 
Remarkably, we find a critical current I c that depends on the 
current flow direction. With more than one transmitting chan- 
nel, the QPC can be tuned to show either a large I c asymmetry 
or a perfect symmetry. In this regime the JJ can act as a su- 
percurrent rectifier jlOlll ill . At the origin of this effect is the 
anomalous supercurrent — proportional to the external field — 
that appears even for zero phase difference <j> between the 
two superconducting leads |2(| 23], Devices based on InAs- 
related materials, that present strong gate-tunable Rashba SO- 
coupling are good candidates to look for these effects [5, 6]. 
In what follows we present the main results of the theory. 
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FIG. 1: (a) Schematic view of the junction. Panels (b) and (c) corre- 
spond to QPC1 and QPC2 simulated with gates placed at z = 30 and 
90 nm on top of the 2DEG, respectively. The conductance with ideal 
2DEG electrodes (dotted line), with the metallic electrodes described 
in the text (thick line) and the current polarization P (thin line) are 
shown. The Rashba coupling strength is a = 20 meVnm. 
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The total Hamiltonian of the system reads 
H = Hqp C + Hp + H L + H c 



(1) 



here Hqpc describes the central 2DEG with the QPC (see 
Fig. 1(a)). In the effective mass approximation Hqpc — 
(pl + P 2 y )/^ m * + a/h(p v a x - p x a y ) + V(x, y) - g[i B S ■ B, 
where the first two terms are the kinetic energy and the Rashba 
SO-coupling, respectively, and V(x, y) is the confinement po- 
tential that defines the QPC. We use a potential that simulates 
the effect of two electrodes held at a distance z from the 2DEG 
02611 . with a gate voltage controlling the height, V g , of the po- 
tential barrier at the center of the QPC. The last term in Hqpc 
is the Zeeman energy. The Hamiltonians Hp and Hp describe 
the right and left superconducting electrodes with an order pa- 
rameter Ap/ L = Aoe ±l ^/ 2 . Finally, He describes the contact 
between the superconductors and the 2DEG. 

For the numerical calculations we discretized the space, 
mapping the Hamiltonian (1) onto a tight-binding-like model. 
We use a square lattice with hopping matrix elements t n and 
ts for the normal (2DEG) and superconducting materials, re- 
spectively lE7ll . The microscopic parameters of the normal re- 
gion correspond to InAs-like materials ||5|]: the effective mass 
is m* — 0.045m e and the electron density n ~ 10 12 cm -2 
(the Fermi energy is Ep ~ 53 meV). We take the total length 
of the junction L = 1.2/im and analyze two point contacts 
denoted as QPC1 and QPC2 corresponding to different val- 
ues of z. We use Ao = 1.5 meV and the coherent length 
£0 = tiVp I Ao = 43 nm corresponding to Nb films, v F is the 
Fermi velocity of the superconductor yl |6[] . 

We evaluate the normal and anomalous propagators that 
contain the information of all the physical quantities of in- 
terest. For Ao = the normal conductance G = J2 a a 1 Go,o> 
is evaluated using the conventional Landauer-like formulation 
026ll . Here G aM i is the contribution to the conductance due to 
incident electrons with spin a that are transmitted with spin 
a'. The spin polarization of the current is defined as P = 
^ (G(,-f —G a i)/G. These quantities characterize the QPC 
in the normal state. For Ao ^ we calculate the Josephson 
current flowing through the right N-S interface, 127 ] 



^)=if E [t l s 3 N^\*M(xi))-t J SN{^(xM{x 3 ))] 

(2) 

here Xi and Xj are coordinates at the edge of the 2DEG and 
the superconducting electrode, respectively, t % g 3 N is the hop- 
ping matrix element connecting neighboring sites at the in- 
terface and the field operator ip^(x) = (ip^ (x) , ip^(x)) creates 

an electron at coordinate x. We choose tg 3 N = (tpj + ts)/2, 
decreasing this value increases the normal scattering at the in- 
terface producing narrow resonances within the central 2DEG 
region with strong influence on the Josephson current 1 1311 . 

The zero-field case: The normal conductance of the system 
presents clear structures on top of the plateaus, shown in Fig. 
1(b) and 1(c) as broad resonances originated in the scattering 
at the electrode-2DEG interfaces. As shown in the same figure 
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FIG. 2: (a) Current-phase relation (CPR) for different gate voltages 
in QPC1. (b) Critical current and conductance as function of V g . The 
points labelled with letters indicate the parameters of the curves of 
panel (a). In (c) and (d) color maps of the y and z magnetization 
calculated with the critical current of curve C of panel (a) are shown. 
The SO strength is a = 20 meV.nm and the lattice parameter is ao = 
3nm. 



these QPCs generate spin-polarized currents with a polariza- 
tion P in the range [0-0.6] depending on the strength of the 
SO-coupling and the number of transmitting channels. 

For superconducting contacts, the current-phase relation 
(CPR) is shown in Fig. 2 for different values of the param- 
eters. We observe two characteristic CPR: resonant-like and 
tunnelling-like (sinusoidal) relations. As V g changes the junc- 
tion alternates between these two behaviors. The critical cur- 
rent I c is defined as the maximum current in the CPR. The 
dependence of I c and the conductance of the normal state G 
on V g is shown in Fig. 2(b). The structure of both curves is 
similar with the peaks located at the same position showing 
that the maximums of I c are due to one-electron resonances 
J3l. The structure of the Andreev spectrum includes a num- 
ber of (dispersive) states with phase-dependent energies and 
(non-dispersive) states confined at each side of the constric- 
tion, details will be presented elsewhere [28]. When there is 
only one transmitting channel, the Josephson current is domi- 
nated by a single (spin) pair of dispersive Andreev states. As 
V g decreases and the QPC opens, the current is the superposi- 
tion of contributions due to different pairs of states. 

Due to the SO-coupling each Andreev state has a well de- 
fined spin texture which in the presence of a supercurrent con- 
tributes to the local magnetization. As an illustration we cal- 
culate the magnetization components (S y (xi)) and (S z (xi)) 
for all the lattice sites within the 2DEG. As shown in Fig. 2(c) 
and 2(d), the supercurrent-induced steady magnetization has 
an in-plane component perpendicular to the current direction 
and an out-of-plane component with the spin-Hall structure 
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FIG. 3: CPR for QPC1 (panels (a) and (c)) and QPC2 (panels (b) 
and (d)) with an external field of g/is-B = 0.3 meV. In (a) and (b) the 
QPC height V g is set to follow a resonance at the first conductance 
plateau for different values of the SO coupling a (in meVnm). In (c) 
and (d) the I((j>) for different V g (in meV) at the second and third 
plateaus are shown. Note the asymmetry between I(<f>) > and 

m < o. 



Effect of external in-plane magnetic fields: The spin tex- 
ture of each Andreev state has a component along the y- 
direction. In general, for a spin-polarizing QPCs, a pair of 
dispersive Andreev states — indicated by |+) and |— ) — have 
(+|S'j,|+) 7^ — (— \S y \— ). Consequently, with an external 
field B in the y-direction the absolute value of their Zeeman 
shifts is different. The effect is illustrated in Fig. 3(a) and 
3(b). In QPCs with a single transmitting channel we follow 
a resonance as the SO-coupling a is increased. For small a 
we observe the characteristic behavior of a resonant state in 
the presence of a field. The distance between the step-like 
changes of the current is a measure of the Zeeman splitting. 
As a increases the Zeeman shift of one of the Andreev states 
decreases and changes sign (the step-like structure crosses the 
<fi = 7r point). In the large a limit (+|S , 2/ |+) « (— \S y \— ) 
and the CPR shows a single step. In this limit, the QPC acts 
as a very efficient spin filter in the first plateau ffl Tl8ll : all 
transmitted electrons have essentially the same spin orienta- 
tion. This results in an anomalous supercurrent for = 0. To 
estimate it, we may assume a smooth QPC and evaluate the 
phase shift i9„ (E) acquired by an electron travelling from one 
superconducting electrode to the other in the WKB approxi- 
mation. In d n (E) the index n = ± is the channel index and E 
is the energy measured from the Fermi energy. For small E 
we have 



w k F X n + - — S n , 



(3) 



here k F (vf) is the Fermi wavevector (velocity) in the 2DEG, 
A n < L and S n > L are related to the effective length of 
the junction at the Fermi energy. Assuming that at resonance 
the scattering at SN interface plays no important role we have 



here ft(x) is a periodic function with fl(x) = x for \x\ < 
7r, A„ = hvF/Sn and p n B is the Zeeman shift. Then, 
in Fig. 3(a) and 3(b) the step-like structures correspond 
to <fi = 7r — p n B/A n and the anomalous current 1(0) = 
(2e/h) J2 n A 4 "^ i s gi ven by the total Zeeman energy of the 
transmitting channels — note it does not depend on L. Non- 
linear effects with larger anomalous currents occur for large 
fields if /j, n B/A n > it for some of the channels. We consider 
only the linear regime in which the current cancels for a phase 
ip = 7r/(0)<5+<5_ /evF(S+ + S-). This (p— junctions in a ring 
geometry generates a spontaneous current with a fraction of 
a vortex threading the ring. Let us point out that the Zeeman 
field couples to the momentum through the SO-coupling, act- 
ing as a gauge field that generate a (^-junction. Yet, this kind 
of symmetry argument is not sufficient, and the adiabatic QPC 
here plays an essential role in filtering and coherently mixing 
very few transverse channels. Indeed, no such effects have 



fJ-n 



B_, 



(4) 



been obtained in wide junctions 12111 and quantum dots 12411 ■ 
An exception is Ref. B23I1 were a ID case was considered. 

As in the zero-field case, a small change in V g shifts the 
resonance from Ep and the current becomes a smooth func- 
tion of the phase, characteristic of a non-resonant junction. 
Remarkably, with more than one pair of transmitting chan- 
nels the CPR presents new effects, displayed for instance in 
the second conductance plateau, Fig. 3(c) and 3(d). In Fig. 
3(c), we show a value of V g for which two resonances — with 
different values of their parameters /i„ and 5 n — lie at Ep- 

Changing V g shifts each resonance by a different amount. 
The total CPR now results from the superposition of contribu- 
tions with different step-like structures and different </?-shifts. 
This leads to a critical current I c that depends on the current 
direction. We define the critical current asymmetry as I* jl~ 
where and I~ are the critical currents for each flow direc- 
tion. Figure 4 shows that for physical values of the param- 
eters a large asymmetry can be obtained by tuning the gate 
voltage. The magnitude of the asymmetry depends on the 
detailed structure of the dispersive Andreev states, which is 
determined by the interplay between the SO-coupling, the ex- 
ternal field and the QPC potential (see figs. 4(c) and 4(d)). We 
found that the asymmetry /+ / 1~ can be larger than 3. These 
large assymmetry values are the main result of our work. 

In summary, we have shown that a spin-polarizing QPC 
brings new physics to the JJs. While the most relevant effect is 
the new mechanism to generate critical current asymmetries, 
the device shows other interesting properties that highlight the 
effects of the SO-coupling on the Andreev states: z) the super- 
current generates a magnetization in the 2DEG, being larger 
at the core of the QPC; ii) an external in-plane magnetic field 
induces an anomalous current at zero phase difference (a in- 
junction) i23ll . In the latter case, we obtain it -junctions for 
some values of the gate voltage (as for a — 0), while in gen- 
eral ip < ir. Such junctions, tunable both by an external flux 
and a Zeeman field, may have applications in SQUIDs or su- 
perconducting quantum bits. With more than one transmit- 
ting channel, the external field induces a large critical current 
asymmetry if the QPC potential V g is properly tuned. Even 
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FIG. 4: Critical current asymmetry vs. V g for QPC1 (a) and QPC2 

(b) and different values of a (in meVnm) with the same applied field 
of Fig. 3. The conductance of the QPCs is also shown (dotted line). 
The lower panels are maps of the supercurrent for QPC 1 with a — 20 

(c) and QPC2 with a — 5 (d) in the [<^,V S ] plane. The V" 9 -scales of 
these maps are shown with horizontal bars in panels (a) and (b). The 
vertical dashed lines correspond to the curves of Fig. 3(c) and 3(d). 



for moderate values of the SO-coupling, and realistic values 
of the external field (B < IT), the asymmetry can be quite 
large. The QPC is a central ingredient as it allows the con- 
trol of the number and properties of the transmitting channels. 
These junctions act as supercurrent rectifiers in the interval 
min(/+,/ ( r) < |/| < ruax(/+, which can be controlled 
by adjusting the gate voltage. As this effect relies on the con- 
trol of the spin polarization of the Andreev states, it gener- 
ates a new alternative for supercurrent rectifiers based on pure 
spintronic effects. 
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